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Abstract 

Here we consider micron-sized samples with any axisymmetric body shape and made with a 
canted antiferromagnet, like hematite or iron borate. We find that its ground state can be a 
magnetic vortex with a topologically non-trivial distribution of the sublattice magnetization / and 
planar coreless vortex-like structure for the net magnetization M. For antiferromagnetic samples in 
the vortex state, in addition to low-frequency modes, we find high-frequency modes with frequencies 
over the range of hundreds of gigahertz, including a mode localized in a region of radius ~ 30-40 
nm near the vortex core. 

PACS numbers: 75.50.Ee, 76.50.+g, 75.30.Ds, 75.10.Hk 
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I. INTRODUCTION 



The magnetic properties of submicron ferromagnetic samples shaped as circular cylinders 
(magnetic dots) are attracting considerable attention mainly due to their potential applica- 
tions (see, e.g., Ref. Circular dots possess an equilibrium magnetic configuration which 
corresponds to a vortex structure just above the single domain length scale, with radius 
R > -Rcrit ~ 100-200 nm. The ferromagnetic vortex state consists of an in-plane flux-closure 
magnetization distribution and a central core with radius ~ 20-30 nm, magnetized perpen- 
dicular to the dot plane. Reducing the magnetostatic energy comes at the cost of a large 
exchange energy near the vortex core, as well as the magnetostatic energy caused by the 
core magnetization. Due to spatial quantization, the magnon modes for dots have a discrete 
spectrum. The possibility to control the localization and interference of magnons spawned 
magnonics, trying to use these modes for developing a new generation of microwave devices 
with submicron active elements (see, e.g., Ref. |2|). 

These samples also provide an ideal experimental system for studying static and especially 
dynamic properties of relatively simple topologically non-trivial magnetic structures which 
are fundamentally interesting objects in different research areas of physics. The investigation 
of the non-uniform states of ordered media with non-trivial topology of the order param- 
eter can be considered as one of the most impressive achievements of modern condensed 
matter physics (see, e.g., Ref. 1^, and field theory (see, e.g., Ref. 4). For example, vortices 
(topological solitons) appear in many systems with continuously degenerate ground states, 
whose properties are determined by some phase-like variable 0, including superconductors 
(see, e.g., Ref. Is]), quantum liquids (hehum-II, and different phases of superfluid ^He, see, 
e.g., Ref. 0), dilute Bose-Einstein condensates,- and also different models of magnets: fer- 
romagnets and antiferromagnets,-""- and spin nematics.— The contribution of topological 
excitations (vortices and vortex pairs) to the thermodynamics and response functions of a 
two-dimensional ordered media are well known. At low temperatures, vortices are bound 
into pairs, forming a Berezinskii phase with absence of long-range order, but with quasi-long- 
range order. The unbinding of the vortex pairs at high enough temperatures, T > Tbkt, 
leads to the Berezinskii-Kosterlitz-Thouless phase transition. In particular, the translational 
motion of vortices leads to a central peak in the dynamic correlation functions, which has 



been observed experimentally; see Refs.191, 
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13 and references therein. 
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Physical systems with a topological defect as a ground state are of special interest. The 
role of vortices in rotating superfiuid systems is known.- Additional examples include single- 
connected samples of the A-phase of superfiuid ^He, where the true energy minimum is a 
non-trivial state (bujoom) with a surface vortex-like singularity.^ Moreover, non-uniform 
states can also appear in small magnetic samples having high enough surface anisotropy.— 
Among these examples, magnets are interesting because samples bearing vortices can be pre- 
pared of micrometer and sub-micrometer size, and vortices can be present at high enough 
temperatures, including room temperature. Furthermore, their static and dynamic proper- 
ties can be observed using different techniques. (e.g., Refs. [l IgI 17) 



A. Summary of results 

All previous studies of magnetic vortices caused by magnetic dipole interactions were 
carried out on magnetic samples made with soft ferromagnets with large magnetization Ms, 
like permalloy, with AnM^ ~ 1 T. In this work, we show that vortices can be a ground state 
for micron-scale samples with axial symmetry, but now made with another important class of 
magnetic materials: antiferromagnets with easy-plane anisotropy and Dzyaloshinskii-Moriya 
interaction (DMI). Some antiferromagnets possess a small, but non-zero, net magnetization 
caused by a weak non-coUinearity of the sublattices (sublattice canting) originated from 
the DMI. As we show here, the formation of a vortex in an AFM sample leads to a more 
effective minimization of the non-local magnetic dipole interaction, compared to the case 
with ferromagnetic dots. Typical antiferromagnets include hematite a-Fe203, iron borate 
FeBOs, and orthoferrites.- These materials exhibit magnetic ordering at high temperatures. 
Moreover, they have unique physical properties: orthoferrites and iron borate are transparent 
in the optical range and have a strong Faraday effect, and the magnetoelastic coupling is 
quite high in hematite and iron borate.^ For antiferromagnets, typical frequencies of spin 
oscillations have values in a wide region, from gigahertz to terahertz, which can be excited 
by different techniques (not only the standard one using of ac-magnetic fields, but with 
optical or acoustic methods as well). Spin oscillations can also be triggered by ultra-short 
laser pulses.—"— For antiferromagnetic (AFM) samples in the vortex state, here we derive 
a rich spectrum of discrete magnon modes with frequencies from sub-gigahertz to hundreds 
of gigahertz, including a mode localized near the vortex core. These results should help to 
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extend the frequency range of magnonic microwave devices till sub-terahertz region. 



II. MODEL 

For antiferromagnets, the exchange interaction between neighboring spins facilitates an 
antiparallel spin orientation, which leads to a structure with two antiparallel magnetic sub- 
lattices, Ml and M2, |Mi| = IM2I = Mq. To describe the structure of this antiferromagnet, 
it is convenient to introduce irreducible combinations of the vectors Mi and M2, the net 
magnetization M = Mi + M2 = 2Mom, and the sublattice magnetization vector / = (Mi — 
M2)/2Mo. The vectors m and / are subject to the constraints: (m ■ /) = and + P = 1. 
As |m| <^ 1, the vector / could be considered as a unit vector. The mutual orientation of the 
sublattices is determined by the sum of the energy of uniform exchange We^ and the DMI 
energy, Wbm, 

We. = H,.Mom^ , VTdm = 2MoHn [d- (m x V)] , 

where the unit vector d is directed along the symmetry axis of the magnet. Here we will not 
discuss the role of external magnetic fields. The parameters iJcx ~ 3-10^-10^ T and Hd ~ 10 
T are the exchange field and DMI field, respectively. Using the energy (Wex + Wbm) and 
the dynamical equations for Mi and M2, one can find,- 

M = Mom x i) + -£ (^i x - j . M„„ = ^ , (1) 

where 7 is the gyromagnetic ratio. The first term gives the static value of the antiferromagnet 
net magnetization Mdm? comprising a small parameter, ifn/-f^ex ~ 10~^, and the second 
term describes the dynamic canting of sublattices; see Ref. l7| for details. Note that Mdm 
is much smaller than either Mq or the value of Mg for typical ferromagnets. However, 
the role of the magnetostatic energy caused by Mdm could be essential, and could lead to 
the appearance of a domain structure for antiferromagnets,-'^^ and the magnetic dipolar 
stabilization of the long-range magnetic order for the 2D case.^^ We will show that for the 
formation of equilibrium vortices, antiferromagnets have some advantages compared with 
soft ferromagnets. 

The static and especially dynamic properties of antiferromagnets are essentially different 
from the ones of ferromagnets. The spin dynamics of an antiferromagnet can be described 
using the so-called sigma-model (a- model), a djTiamical equation only for the vector / (see. 



e.g., Ref. [7). In this approach, the magnetization M = Mi + M2 = 2Mom is a slave variable 
and can be written in terms of the vector / and its time derivative, see Eq. ([1]). Within the 
cr-model, the equation for the normalized (unit) antiferromagnetic vector / = L/\L\ can be 
written through the variation of the Lagrangian C[l] 
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2c2./ \dt 



d\-W 



(2) 



where A is the non-uniform exchange constant, c = AH^y^IMq is the characteristic speed 
describing the AFM spin dynamics, [/] is the functional describing the static energy of the 
AFM. It is convenient to present the energy functional in the form = PVo[^ + W^m, 



MH„,dx- 



(3) 



where the first term PVo[^] determines the local model of an AFM, including the energy 
of non-uniform exchange and the easy-plane anisotropy energy through only the vector /, 
> is the anisotropy constant, and the xy— plane is the easy plane for spins. The second 
(non-local) term Wm is the magnetic dipole energy, Hm is the demagnetization field caused 
by the AFM magnetization ([1]), the field Hm is determined by the magnetostatic equation, 

div(#„ + 47rM) = 0, curl^^ = , 

— * 

M = 2Mom, with the standard boundary conditions: the continuity of the normal compo- 
nent of {Hm + 4:71 M) and the tangential component of Hm, on the border of the sample. 
Thus, sources of Hm can be considered as formal "magnetic charges" , where volume charges 
equal to divM and surface charges equal —M ■ n, where n is the unit vector normal to the 



border ( see the monographs 
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for application 



for general considerations and Refs. 
to vortices). As mentioned above, the presence of surface anisotropy with the constant -ft'surf 
can contribute to the energy of a micron-scale sample leading to non-uniform states.— 
But this contribution is essential for a high enough anisotropy Ksuri ^ K and quite small 
samples, and we will neglect this contribution, and concentrate on an alternative source of 
non-uniform states, namely, the magnetic dipole interaction. 
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A. Dynamic sigma-model 



The variation of the total Lagrangian ([2]) gives a dynamic cj-model equation, where the 
static spin distribution for an uniaxial antiferromagnet is determined by the minimization 
of an energy functional of the form of Eq. It is useful to start with a variation of the 
local energy Wq[1] only, that gives a general two-dimensional (2D) vortex solution for the 
vectors Z, M of the form 

I = Cz cos9 + sm9[exCos{x + fo) + (^y^^^ix + fo)] y 
M = Mdm • sin 9[ey cos(x + V^o) - 4 sin(x + v^o)], (4) 

where 6 = 6{r), r and x ^I's polar coordinates in an easy-plane of the magnet, the vector 
is the hard axis, and the value of ipo is arbitrary The function 6{r) is determined by the 
ordinary differential equation 

'^de . , , f 1 \ r 
+ -X: = sin^cos^ — -1 ,x = -, (5) 



dx'^ X dx \ / ' Iq 

and exponentially tends to 7r/2 for r ^ Iq, with characteristic size Iq = \/ AjK. Also, in the 
center of the vortex (at r = 0), the value of sin 6'(r = 0) = 0, see Fig. 1. 

In the region of the vortex core, / deviates from the easy-plane, and the anisotropy 
energy increases. The state (jl]) is non-uniform, increasing the exchange energy. Therefore, 
for the local easy-plane model, the appearance of a vortex costs some energy, i.e., the vortex 
corresponds to excited states of antiferromagnets. Also, vortex excitations are important 
for describing of thermodynamics of 2D antiferromagnets.— As we will show in the next 
section, considering the dipolar energy a vortex can be the ground state of a circular 

magnetic sample made with canted antiferromagnets. 



III. ENERGY OF THE VORTEX AND UNIFORM STATES 

For small samples made with canted antiferromagnets, the energy loss caused by a vortex 
within the local model W = Wo[l] can be compensated by the energy of the magnetic dipole 
interaction. To explain this, note that for a uniform distribution state, the contribution Wm 
unavoidably results in a loss of the system energy, which is proportional to the volume V of 
the sample.—"^ The energy of the uniform state could be estimated as 

^(homog) _ 2nNMl^V = 2TiNMl{2H^/H^fV, (6) 
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FIG. 1: (Color online) Schematic structure of the AFM vortex for cos 9{r = 0) = 1 and the special 
value of the arbitrary parameter ipo = 0. The vectors / (thin blue arrows in the radial direction) 
and M (thick tangential red arrows; shown not in scale) are depicted for the area of the core 
(dashed circle) and far from it ( larger dotted circle). The red dot at the origin indicates the value 
M = for the state with / = Sz perpendicular to the plane. 

where N is the effective demagnetizing factor in the direction perpendicular to the sample 

22 24 

axis.— 

For magnetic samples other than an ellipsoid (for example, for a cylinder with a finite as- 
pect ratio R/L), the distribution of the field is non- uniform, leading to a non-topological 
quasi-uniform state known as a flower state.— But a detailed numerical analysis shows,—"— 
in all the stability region of the quasi-uniform state where Eq. (jG]) is valid, the correspond- 
ing value of N is practically the same as for the uniform magnetization, and the numerical 
results^E"^ obtained within this approximation are in good agreement with experiments.— 

For the static vortex state (j4]) with a chosen value of sin<^o = 0, with M oc {d x I), one 
can find 

M = ctMdm ■ sin6'(-ea; sinx + cosx), (7) 

where a = cosv?o = ±1- A unique property of the state ([7]) is that it can also exactly 
minimize the energy of the magnetic dipole interaction Wm, giving Hm = in the overall 
space. Indeed, the projection of M on the lateral surface of any axisymmetric body with its 
symmetry axis parallel to the ^-axis equals to zero. Also divM = for the distribution given 
by Eq. ([7]). Moreover, due to the symmetry of the DMI, (M ■ d) = 0, the distribution of the 
magnetization M (JTj) is purely planar (in contrast to l) and the out-of-plane component of 
M is absent. In the vicinity of the vortex core, the length of the vector M decreases, turning 
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to zero in the vortex center, see Fig. 1. Such feature is well known for domain walls in some 
orthoferrites.- Thus, the AFM vortex is the unique spin configuration which does not create 
a demagnetization field in a singly-connected body (for ferromagnet with \M\ = const, a 
configuration with = can be possible only for magnetic rings having the topology of a 
torus) . 



A. Comparing the energies of the vortex and uniform states for antiferromagnets 

Let us now compare the energies of the vortex state and the uniform state for AFM 
samples shaped as a cylinder with height L and radius R. For the vortex state, Hm = 
over the the sample volume, and the vortex energy is determined by the simple formula,- 

E, = 7rAL\n(p^], (8) 



\ ^0/ 

where p 4.1 is a numerical factor. For long cylinders with L ^ R, the value of ^ 1/2 
and the vortex state becomes favorable if the radius R exceeds some critical value -Rent, 



R > Rrr\t = 2/dip4 / In ( -T-^ 







where /dip determines the spatial scale corresponding to the magnetic dipole interaction. 
Note that Zdip comprises a large parameter He/ Hd ~ 30-100, and /dip ^ h- In the case of a 
thin disk, L <^ R and the demagnetization field energy becomes^i^^ 

^(homog) _ 27ri?L'M^^ln(4i?/L), 

and the vortex state is energetically favorable for RL > (RL) cnt = 2/^ip. 

For concrete estimates we take the parameters of iron borate, A = 0.7 • 10^® erg/cm, 
A' = 4.9 ■ 10*^ erg/cm^ and AttMbm = 120 Oe. Then we obtain that Iq = 3.8 nm, i.e., 
the core size is of the same order of magnitude as for a typical ferromagnet (for permalloy 
Iq = 4.8 nm). The value /dip is essentially higher, e.g., for iron borate /dip = 220 nm. 
Combining these data one finds for the long cylinder Rent = 0.9 fim. For a thin disk sample 
the characteristic scale has submicron value: ^y (RL) cnt = 0.4 /im. Similar estimates are 
obtained for orthoferrites, and somewhat higher values for hematite. Thus, despite the fact 
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that the characteristic values for the dipole length /dip for a ferromagnet and antiferromagnet 
differ hundredfold, the characteristic critical sizes differ not so drastically (for permalloy 
-Rcrit ~ 100-200 nm). This is caused by the aforementioned fact that the magnetic field 
created by the vortex core is completely absent for the AFM vortex. The situation here is 
common to that for ferromagnetic nanorings, where the vortex core is absent. Despite the 
fact that the vortex core size in ferromagnetic dots is rather small, the core contribution to 
Wm for ferromagnetic dots of rather big radius R > 0.5 fim is negligible, but it becomes 
essential for small samples with R close to the critical size. Note as well that the vortex core 
magnetic field in the ferromagnet destroys the purely 2D character of the distribution of M 
in Eq. (jlj), and the core size changes over the thickness of the sample. For an AFM vortex, 
the value of equals exactly zero, and a truly 2D distribution of I and M, independent 
of the coordinate z along the body axis, is possible. 

IV. MAGNON MODES FOR VORTEX STATE ANTIFERROMAGNETIC SAM- 
PLES 

The dynamics of / considered within the a-model approach differs from the dynamics of 
a ferromagnetic magnetization described by the Landau-Lifshitz equation. The main dif- 
ference is that the cr-model equation contains a dynamical term with a second-order time 
derivative of /, whereas the Landau-Lifshitz equation is first-order in time. For this reason, 
for antiferromagnets, two magnon branches exist, instead of one for ferromagnets<- For both 
AFM modes, the elliptic polarization of the oscillations of Mi and M2 is such that the oscil- 
lations of the vector / have a linear polarization.-^ For an easy-plane antiferromagnet, these 
two branches are a low-frequency quasi-ferromagnetic (QFM) branch and a high-frequency 
quasi-antiferromagnetic (QAF) branch, respectively. QFM magnons involve oscillations of 
the vectors I and M in the easy-plane, with a weak deviation of M from the easy-plane 
caused by the last term in ([1]). The second QAF branch corresponds to the out-of-plane 
oscillations of I with the dispersion law 



where k is the magnon wave vector. The gap of the QAF branch, Ug, contains a large 
value Hex and attains hundreds of GHz. Thus compared to ferromagnets, both the magnon 
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frequency and speed for AFM dynamics contain a large parameter 

^HJH, ~ 30-100, 

(iJex and iJa are the exchange field and the anisotropy field, respectively) which can be 
referred as the exchange amplification of the dynamical parameters of AFM. The frequency 
Wg of AFM magnon modes reaches hundreds of GHz, with values ~170 GHz for hematite, 
100-500 GHz for different orthoferrites, and 310 GHz for iron borate.— Recent studies using 
ultra-short laser pulses showed the possibility to excite spin oscillations of non-small ampli- 
tude for orthoferrites^ii^ and iron borate.-^ This technique can be also extended to other 
antiferromagnets, including those without Dzyaloshinskii-Moriya interaction.— 

Since the magnon spectra of bulk ferromagnets and antiferromagnets differ significantly, 
one can expect an essential difference for the magnon modes of the vortex state for both 
AFM and ferromagnetic samples. Let us briefly recall the properties of normal modes for 
disk-shaped vortex state ferromagnetic samples (ferromagnetic dots). For such dots, the 
presence of a discrete spectrum of magnon modes, characterized by the principal number 
(the number of nodes) n and the azimuthal number m, is well established.—"—'' This spectrum 
includes a single low-frequency mode of precessional motion of a vortex core (n = 0, m = 1) 
with a resonant frequency in the subGHz region,— a set of radially symmetrical modes 
with m = 0,"^^ and also a system of slightly splinted doublets with azimuthal numbers 
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m = ±|m|, with frequencies u\m\,n 7^ <^-\m\,n, and uj\m\,n — ^~\m\,n ^ ^\m\,n] see Refs 
The same classification is valid for vortices for local easy-plane ferromagnets.— Wysin^ 
demonstrated the direct correspondence of the gyroscopic character of vortex dynamics and 
doublet splinting. 

For an AFM in the vortex state, small sample each of two magnon branches, QFM and 
QAF, produce a set of discrete modes with given n and m; however, their properties are 
different compared to that of a ferromagnetic dot. Below we present a general analysis of 
small oscillations above the vortex ground state. 

A. General equations and mode symmetry. 

The dynamics of small deviations from the AFM static vortex solution will be considered 
here for a thin circular sample (AFM dot) only, where the z-dependence of the vectors / 
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and m can be neglected. It is convenient to introduce a local set of orthogonal unit vectors 
6*1, 6*2 and 6*3, where 63 coincides with the local direction of the unit vector / in the vortex, 
63 = l{x, y) = cos 6*06^ + sin 9or, see Fig. [H 62 = —e^ sin x + cos x = X; and Ci = (e2 x 63). 
It is easy to see that the projection (/ ■ ei) = i} describes small deviations of 6 from the 
vortex ground state, 9 = Oq (r) + {}, where 9q is the solution of Eq. (E]), describing a static 
vortex structure, and fi = {I ■ 6*2) is the azimuthal component of /. In linear approximation, 
the equations for ?9 and /i become a set of coupled partial differential equations, 

2cos^o<9/i _ /q^^^ 



X 



2 dx dt'^ ' 



^2 T./ N r T 1 2cos6'o(9'(9 

v^-ViW-ft]^-— = (10) 

where x = r /Iq, 9o = 9o{x) is the solution of Eq. (j5]), and = ZqV. The equations are 
symmetric in ^ and yU with local Schrodinger-type differential operators in front, as well as 
with non-local parts Hi, H2- The "potentials" Vi, V2 in local Schrodinger-type operators are 
determined by Wq [/] , see Eq. (|3]) ; they have the same form as for easy-plane magnets,- 

1 A 2 a fd9o^ ^ 



_41j43 



V,= ^-^-ljoos2»„.l/.= ^^-ljcos^»„-^-^j . (11) 

Note that the potentials in this Schrodinger-type operators are not small, but localized near 
the vortex core. Non-local magnetostatic effects, defined by magnetic dipole interactions, 
are included in the integral operators. Hi, H2- For a ferromagnetic vortex, their form was 



determined and their role was discussed in Refs. I26U36I . Generally, for antiferromagnets 

— * 

the magnetization includes not only terms proportional to the in-plane components of I, 
but also time derivatives of the vector /. For this reason, the structure of these operators 
presented though the vector / is much more complicated than the corresponding structure for 
ferromagnets. But the non-local contributions to f lTOj) are essential only for in-plane modes 
with low frequencies uj <^ 'jHdm, and for this case the dynamical part of the magnetization 
is negligible, as shown below. By means of this approximation, one can demonstrate that the 
operators Hi^2, corresponding to volume and to edge magnetic charges, take the same form 
as for a ferromagnetic vortex, after replacing Mg — )■ Mom and — )■ 0. In particular, the 
angular dependence of the eigenfunctions for these operators is the same as for ferromagnetic 
vortices, namely, i^i^2 exp(zmx) = exp(imx), where the integer m (m = 0, ±1 ±2, . . .) 
is the azimuthal number. Thus for AFM vortices, even considering the non-local magnetic 
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dipole interaction with M oc [l x d), the separation of the radial and azimuthal parts of the 
deviations is possible, and the magnon modes are of the form exp{imx) fm{r) ■ This property 
is of importance for the problem of magnon modes above the AFM vortex ground state. 

Thus, the static part of the equations ( ITOl) . both local and non-local terms, have the same 
form as for the well-studied case of the ferromagnetic vortex, but the dynamical parts differ 
strongly. This produces a crucial difference in the magnon modes of these magnets. For 
ferromagnetic vortices, the magnon eigenstates {{^m , f^m} depend on x and t in combinations 
as sin(mx + (^t) or cos(mx -|- ut), whereas for a AFM vortex, a more general ansatz of the 
form 

= /a(r)(Ae*™^ + Be-'"''') exp(zw„t) + c.c. 
/ia = i^a(r)(Ae*™^-5e-^™^)exp(ia;„t) + c.c. (12) 

is appropriate.— Here a = (n, m) is a full set of discrete numbers labelling the magnon 
eigenstates, and n is the nodal number. Substitution of this ansatz demonstrates, in contrast 
to the case of a ferromagnetic vortex, the full degeneracy of the frequency over the sign of 
m.— As a change of sign in the number m can also be interpreted as a change of the sense 
of rotation of the eigenmode (change of sign in the eigenfrequency u), physically we have 
the situation of two independent oscillators rotating clock- and counter-clockwise with the 
same frequency (which can also be combined to give two linear oscillators in independent 
directions). This degeneracy was clearly demonstrated by solving the ordinary differential 
equations for fa and ga ( ITOl) . as well as due to direct numerical simulations of the magnon 
modes above an AFM vortex.^^ Thus, the absence of gyroscopical properties for the a-model 
equation is manifested in the fact that for a AFM vortex the modes with azimuthal numbers 
m = \m\ m = —\m\ are degenerate, i.e., the splitting of doublets with m = ±|m|, typical 
for the ferromagnetic vortex, is completely absent.— 

Note one more important difference from the ferromagnetic case: for the AFM vor- 
tex the coupling of in-plane and out-of-plane oscillations comes only from the term with 
{cos6o){d/dx)- This means, that (i) for any mode the coupling vanishes exponentially at 
X — 7- oo; (ii) the modes of radially symmetric (m = 0) in-plane and out-of-plane oscillations 
are completely uncoupled. Both properties will be employed below for calculating of the 
magnon frequencies. 
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B. Collective variables for vortex core oscillations. 



First note that the equation ( ITOj) for an infinite magnet has a simple zero-frequency 
solution with m = 1 

-d = {a ■ V6'o) , = sin6'o(a ■ V0o)) 

where a is an arbitrary vector, and 60, (po describe the static vortex solution. Indeed, this 
perturbation describes a displacement of the vortex for a (small) vector a. Such "zero 
mode" appears for any soliton problem, refiecting an arbitrary choice of the soliton (vortex) 
position.- For finite-size magnets, such modes beget low-frequency modes corresponding to 
the motion of a vortex core. Their analysis with the equations ( ITOl) is quite complicated, but 
it can be done within the approach based on the scattering amplitude formalism, which is 
developed for the Gross-Pitaevski equation^ and local models for easy plane magnets.^- But 
there is an easier and more convenient way to calculate the frequency of this mode based 
on a collective variable approach. Here the collective-variable is the vortex coordinate X, 
which motion is described by a characteristic dynamic equation.— 1^ 

Thus, for an AFM vortex, as well as for ferromagnets, one can expect the appearance 
of a special mode of vortex core oscillations. But the dynamical equations for the AFM 
vortex core coordinate differ significantly from that for ferromagnets. The a-model equation 
contains a dynamical term with a second time derivative of /, combined with gradients of 
/ in the Lorentz-invariant form (Pl/dt^ — c^V^/, whereas the Landau-Lifshitz equation is 
first-order in time. The chosen speed c = '-)^/AI^JMq plays roles of both the magnon 
speed and the speed limit of solitons, it is only determined by the exchange interaction and 
attains tens km/s; e.g., c ~ 1.4-10^ m/s for iron borate and c ~ 2 - 10^ m/s for orthoferrites.- 

The formal Lorentz-invariance of spin dynamics of antiferromagnets manifests itself in 
the motion of any AFM solitons, in particular, the motion of the AFM vortex core:— the 
dynamical equation for the core coordinate at small vortex speed X, when at \dX/dt\ <^ c, 
possess an inertial term, 

where the effective vortex mass Mv = ^v/c^ and F is an external force acting on vortex. 
For the case of interest here, the free dynamics of the vortex in a circular sample, F is 
the restoring force: in linear approximation F = —nX, where k is the stiffness coefficient. 
With this force, the vortex core dynamics is not a precession, as for the gyroscopic Thiele 
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equation for ferromagnetic vortices,—"— but rectilinear oscillations, X(t) = acos^u^t + (p^), 
degenerate with respect to the direction a and 0o, with frequency lo^ = a/ k/M^. For the 
easy-plane AFM model with Wm = 0, such dynamics has been observed by direct numerical 
simulations.'^'^ For a vortex state dot with R > Rent, the value of k is determined by the 
demagnetizing field, its value can be obtained from the known value for a ferromagnet by 
replacing Mg Mdm,— which gives k = 10 ■ AnM^^L'^ /9R, and 

2cM^mVT0L 
3^^AR\ii{pR/lo) ■ 

A simple estimate gives that u^, as for a ferromagnetic vortex, is in the subGHz region, 
but with different (approximately square root, instead of linear for ferromagnetic vortex) 
dependence on the aspect ratio L/R. 



C. Other low-frequency modes. 

Far from the vortex core, the other modes from this set are approximately characterized 
by in-plane oscillations of / and M. As their frequencies are small, u ^ 'jHdu, for these 
modes the magnetization M is determined mainly by the in-plane static contribution ([1]), 
and the formulae for the demagnetization field energy for ferromagnetic vortices can be 
used. Moreover, the data known for the ferromagnetic magnon frequencies u™ (in first 
approximation over the small parameter L/R, i.e., in the magnetostatic approximation), can 
be directly used for the calculation of the corresponding frequencies for magnon modes above 
the AFM vortex ground state. Note that, in this approximation, = = u™, because 
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the doublet splitting for a ferromagnetic vortex is proportional to {L/RY; see Refs. 

To make this essential simplification, note that for the ferromagnetic case the Landau 
Lifshitz equation, linearized over the vortex ground state, can be written as dm^/dt = 
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Here the dimen- 



Air'-fMsTJiz, drriz/dt = h2Msmr, as was shown in Appendix B of Ref. 
sionless operator h2 determines the non-local magnetostatic part of the magnetic dipole 
interaction. Then the equations can be easily rewritten as d'^rrir/dt^ + (ci;™)^mr. = 0, 



where = Att'^Ms\J {h2), where (/i2) is the eigenvalue of the operator h2. These values 
for modes with different angular dependence can be either estimated theoretically or taken 
from experiments.— 

Now we will return to the case of AFM vortices. Neglecting the vortex core contribution, 
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the second equation of the system (fTOj) reduces to the form 

aV'9t' = (c//o)'i^2/i, 

having exactly the same structure as the equation for for a ferromagnetic vortex. For 
the simphfied form of the magnetization, M = Mdm(^ x Cz) the same "magnetic charges", 
both volume charges, divM, and surface charges, (e^ ■ M), as for a ferromagnetic vortex, are 
produced. Thus the operator h2 differs from H2 by a simple scaling relation, Mg — ?• Mdm and 
A/AttM^ ll, which gives H2 = {AixMl^j/ K)h2. Then, it is easy to obtain the frequency 
of the AFM mode Um in terms of the frequency of the ferromagnetic mode a;™ (if it is 
known) as follows 

'""^ - 7M.y4^ ■ ^^^^ 
In particular, the frequency of radially symmetric oscillations, having the highest fre- 
quency — for modes with minimal nodal number n, can be presented through the known 
value for as^ 



= Jl„ . (15) 



For these modes, the frequencies are of the order of a few of GHz, with an approximately 
square root dependence on the aspect ratio L/R. 

Note the absence of gyroscopic properties for the a-model equation. This is manifested 
both in the absence of a gyroforce for an AFM vortex, as well as in the fact that for an AFM 
vortex the modes with the azimuthal numbers m = \m\ m = — |m| are degenerate. The 
splitting of doublets with m = ±|m|, typical for ferromagnetic vortices, is absent for AFM 
vortices.—"^ 



D. Out of plane high-frequency modes. 

For an AFM sample in the vortex state, the high-frequency QAF branch of magnons 
begets a set of discrete modes with frequencies of the order of Ug, i.e. hundreds of GHz. For 
all these modes far from the vortex core, oscillations of the vector / are out of plane, and the 
static contribution to the weak magnetization, mgtatic oc (e^ x /) is absent, see ([1]). Moreover, 
it is easy to show that the dynamical part of m, mdyn oc {Ixdl/ dt) has a vortex-like structure 
and does not lead to magnetic poles neither on the up and down surfaces nor at the edge. 
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Hence, the dipole interaction is not essential for the description of high-frequency magnons, 
and the results obtained earlier for the vortex in easy-plane antiferromagnets without the 
magnetic dipole interaction—'^^ can be used as a good approximation. In particular, the 
frequencies ujn,m are close to Ug, and the difference {un,m — ^g) decreases as the dot radius 
increases. 

For an AFM vortex within the local easy-plane model, the set of high-frequency radially- 
symmetric modes with m = includes a truly local mode with an amplitude exponentially 
localized within an area of the order of 5/o and with a frequency ui ~ 0.95u;g independent 
on R. Note that the frequency of this mode is inside the range of low-frequency in-plane 
modes, and hence, we have an example of a truly local mode inside a continuum spectrum. 

The presence of a local mode inside the frequency region of low-frequency modes is quite 
a delicate feature, and it is interesting to discuss whether or not such mode survives for the 
vortex state AFM dots when accounting for dipole interactions. It is easy to show that for 
these modes the oscillations of / have no in-plane component, only the ^ ^ 0, even inside 
the core region. Thus, for this mode the oscillations of fh are such that they do not disturb 
the vortex-like closed-flux structure 

6m oc e^[{d'd/dt) + ^H-^'d cosOq] , 

and all the magnetic poles vanish exactly. Hence, for cylindrical dots made with canted 
antiferromagnets in the vortex state, a radially symmetric mode with exponential localization 
inside an area of radius 30-40 nm near the vortex core appears. The frequency of this mode 
is approximately 5% below the energy gap of out-of-plane modes Wg, which gives ~ 9 GHz 
for hematite and ~ 15 GHz for iron borate. We now stress that such modes are absent for 
ferromagnetic vortices. This mode can be imaged as an oscillation of the vortex core size, 
with keeping the in-plane vortex-like structure for m. The total magnetic moment connected 
to this oscillations is zero, and the excitation of these oscillations by a uniform magnetic 
field, either pulsed or periodic, is impossible. However, such oscillations can be excited by an 
instant change of the uniaxial anisotropy, which determines the vortex core size. The novel 
technique^"— of spin excitations by ultra-short laser pulses can be applied here, because the 
linearly polarized light at inclined incidence, due to an ultra-fast inverse Cotton-Mouton or 
inverse Voigt effect, is equivalent to the necessary change of uniaxial anisotropy. 
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V. CONCLUSION. 



To conclude, for micron-sized samples of typical canted antiferromagnets, their ground 
state exhibits a topologically non-trivial spin distribution. The magnetizations of each sub- 
lattice Ml and M2 are characterized by a vortex state with a standard out-of-plane structure, 
but the net magnetization M = Mi + M2 forms a planar vortex, where the projection of the 
magnetization normal to the vortex plane everywhere in the sample is zero; in particular, 
M = in the vortex center. The vortex state AFM dots possess a rich variety of normal 
magnon modes, from rectilinear oscillations of the vortex core position with sub-GHz fre- 
quency to out-of-plane modes with frequencies of the order of hundreds of GHz, including 
a truly local mode. The use of QAF modes for vortex state AFM dots, particularly the 
truly local mode, would allow the application of magnonics for higher frequencies until ~ 
0.3 THz. This mode can be excited by ultra-short laser pulses with linearly polarized light. 

Our theory would be applicable to other systems with an AFM spin structure, like a 
ferromagnetic bilayer dot containing two thin ferromagnetic films with an AFM interaction 
between them, described by the field ifcx- If -f^cx is large enough, ifcx > 47rMs, the anti- 
phase oscillations of the magnetic moments of the layers produce high-frequency modes with 
frequencies of the order of ^7i7cxi^™, where are the frequencies of the modes for a 
single layer dot. 
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